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Fourteen years ago it was suggested that driven diffusive systems may provide suitable 
models for fast ionic conductors and solid electrolytes [[[]. This triggered the interest in 
such systems, which over the years have become a very good scenario to investigate the 
properties of far from equilibrium steady states 0. A simple microscopic realization of 
a driven diffusive system is an Ising lattice gas in contact with a thermal bath in which 
an external uniform driving force is introduced Q. The usual hopping rate to nearest 
neighbours is modified through a bias in the rates for hopping along the direction of the 
external field. In a finite lattice with open boundary conditions, the gas evolves to a non- 
equilibrium steady state with a non-zero particle current. Due to the energy exchange of 
the system with the heat bath and the driving field, the system also has a steady energy 
current. So far driven diffusion of scalar quantities, like the excess of transported particles, 
has been considered. The purpose of this paper is to extend these considerations to vector 
quantities. 

We begin by very briefly recalling the formulation for the standard case of a scalar 
quantity 0. Space is taken to be d- dimensional and parameterized by coordinates (xi) = 
(x||,x_l), with || the direction of the driving force and _L the directions orthogonal to it. 
The deviation s(x, t) of the density from its uniform average satisfies a continuity equation 
d t s + dtJi(s) = T], with Ji(s) the current generated by ordinary diffusion and the drift 
in the direction of the driving force, and r)(x, t) the noise accounting for the heat bath. 
The contribution of ordinary diffusion to the current is proportional to the gradient of s, 
while the drift in the || -direction can be estimated, for small density fluctuations s, by 
expanding in powers of s and retaining terms through order two. This gives the driven 
diffusion equation || 




In the first papers the noise was assumed to be the derivative of a random current, thus 
conserving the number of particles. In subsequent papers [f§ gaussian noise (describing 
systems with random particle sources and drains and preserving particle number only in 
the mean) was also considered. For both types of noise it was shown that driven diffusive 
systems have, for dimensions less than the critical dimension, infrared (IR) stable fixed 
points at which various quantities scale. Furthermore the corresponding critical exponents 
have been computed using the e-expansion || [| |, the values obtained showing that 
density fluctuations spread faster than diffusively. 

It is worth noting that eq. flU) also describes the toppling of grains of sand in a sandpile 
0. As a point of fact, the system described by eq. ([!]) with gaussian noise was first 
considered as a model for the hydrodynamic regime of a running sandpile. However one 
must mention that, after a shift, gaussian noise has vanishing average and this requires 
particles to be randomly added and removed, whereas in the running sandpile particles 
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are only added. In other words, gaussian noise does not reproduce the feeding rate of the 
running sandpile ||. It is also interesting to mention that for d — 1 a field redefinition 
maps the driven diffusion equation (Q) into the Kardar-Parisi-Zhang equation [[|, which 
in turn is identical to the Navier-Stokes equation. This is no longer true for d ^ 1, for 
which the three equations are different. 

We now generalize driven diffusive systems to the case where the external forces acting 
on the system produce drift in a hyperplane of dimension n parameterized by coordinates 
x a (a = 1, . . . , n). We do this by increasing the number of field variables from 1 to n, 
so that we have a vector s a (x,t). The component s a (x,t) describes the excess of grains 
transported along the direction x a . For every component s a (x,t) we write a continuity 
equation d t s a + d{Ji a = r\ a . Proceeding now as for the the standard case, we obtain 

d t s a - ad p dps a - 2cr£ d a dpsp - u±A ± s a + ^ d a s 2 + k dp (ns a s p - 5 a ps 2 ) = r] a , (2) 

where we have used the notation s 2 = s a s a . This equation has two anisotropy coefficients, 
a and <t£, and two coupling constants, A and k. For later convenience, we introduce the 
coefficient 

u = a(l + 20 . (3) 

Eq. (0) makes sense for 1 < n < d. If n — 1, we recover the standard driven diffusion 
equation ([!]); if n — d, anisotropy is maximal and z/j_ = 0. We want to investigate the 
dynamical critical properties of driven diffusive systems based on this equation. In this 
paper, we concentrate on gaussian noise 

( Va {x,t))=0 

(4) 

( Va (x, t) rjpix', t')) = 2 5 af3 ( x - x') 5{t - tf) 

and use the renormalization group and the e-expansion to characterize to first order in e 
the universality classes of the model for dimensions below the critical dimension. A similar 
program can be carried for conserving noise. 

We have already mentioned that eq. (0) may be used to describe drift in a hyperplane. 
Scenarios in which this may become relevant are oil diffusion and diffusion in disorder 
media. More generally, eq. (§) may be considered in connection with directed percolation 
0. We recall that in directed percolation, when the probability p of jumping to the nearest 
neighbour along a preferred direction is set equal to a critical value p*, a phase transition 
occurs. At this transition, the system shows scaling behaviour together with non-standard 
features, like e.g. fractal dimensionality for the percolating cluster. One could think of 
the parameter £ as playing the role of the probability p, with the difference that £ is 
not a priori set equal to any critical value but enters the model as one other parameter. 
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The continuity equation (|J) can also be thought of as describing the toppling of grains 
of sand in an n-dimensional sandpile, with s a (x,t) the deviation in the direction x a of 
the sandpile profile from its uniform average. Since in general n / 1, one may think of 
using eq. (fj) to describe the hydrodynamic regime of sandpiles with avalanches advancing 
anisotropically. Yet we must keep in mind that gaussian noise does not describe the feeding 
rate of a sandpile. Other Langevin equations based on gaussian noise describing surface 
sandpiles have been proposed in ref. [ 10[| . They differ from @ in that they double the 



number of field variables but keep the number of transport directions equal to one. It is 
important to realize that eq. (0) is not the spherical generalization of the ordinary driven 
diffusion equation ([3]) in which, while keeping the number of anisotropy directions equal 
to one, s(x,t) becomes an n- vector field (n = 2j + l) carrying an n-dimensional irreducible 
unitary representation of 0(3). For isotropic spherical generalizations of the Navier-Stokes 
equation and the Kardar-Parisi-Zhan equation, see refs. M and M. 



To study the system (@) + (f|), we use ref. [11] and recast it as a stochastic quantum 
field theory model with classical action 



S„ = dtdx 7T n 



- TT a + d t s a - adpdps a - 2cx£ d a dpsp - z/ ± Aj_s c 
+ ^ d a s 2 + Kdf3 (ns a S{3 - 5 af3 s 2 ) 



(5) 



where 7i a (x, t) is the response field [Q to the assumed gaussian noise. The action S n has 
the following symmetries: 

(a) Translational invariance in all coordinates, 

(b) 0(n) ® 0(d—n) invariance with respect to spatial coordinates, 

(c) Joint reflection of longitudinal coordinates, field variables and response fields, 

(d) Scaling of longitudinal coordinates (t,x a ,^±) — > (t, gx a , xj_), with the anisotropy 
coefficients, coupling constants and fields transforming as (a, £, v±) — > (g 2 ^,^,^), 
(A, k) — > £( n + 2 )/ 2 (A, k) and (s a , ir a ) —> g~ n / 2 (s Q , ir a ). Since the parameter £ is scale 
invariant, we introduce the invariant coupling constants 

1 A Ik 

u = , , „, ,„ v = 



47T (j( n + 2 )/ 4 ~ 47T (j( n + 2 )/ 4 ' 

In addition, if n = 1, the action is invariant under Galilei transformations [|j] 

s(x, t) — > s(xj_, x\\ — Xat, t) + a 
7i(x,t) — > 7r(xj_,a;|| — \at,t) . 

We note that this symmetry does not generalize to n > 1. The reason is that Galilei's 
invariance requires a constant velocity vector and this is only possible for n = 1, since in 
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this case there is only one preferred direction. In other words, it is precisely the passage 
from n = 1 to n ^ 1 what destroys the Galilei symmetry. 

The critical dimension of the model is d c = 4. For d = 4 the model is ultraviolet 
renormalizable by power counting, the only primitively divergent one-particle irreducible 
Green functions being the self-energy (vr Q (p)s ( g(— p)) and the vertex (vr a (p)s ( g(g)s 7 (— p — q)), 
which have linear and logarithmic degrees of divergence respectively. Using dimensional 
regularization |13|], we obtain the following one-loop renormalization group equations: 

du 



Pu 
Pv 

Pi 
c 



d\nfi 

dv 
dln.fi 

d£ 
dln.fi 

din a 
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. e „, n + 2 r 
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(6) 
(7) 
(8) 
(9) 
(10) 



where fi is the dimensional regularization mass scale, e = 4 — d is the deviation from the 
critical dimension and J 7 , Q, Ti and £ are functions of u, v and £ given by 

n (n + 2) X(u, v, £) = u 2 x uu (£) + uv x uv (£) + v 2 x m {£) X = J 7 , Q, H, £ . 

For completeness we write the expressions of x uu , x uv and x vv \ 
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with a±, 02, a 3 and a 4 taking the form 
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and -F(a, b; c; 2) being the hypergeometric function. For n — 1 there is no K-term in so 
that k = and hence v = 0. From this, and using n = 1, it follows that J 7 = 7i + nC = 0. 
Thus f3 u = e/2 and, as a result, /3 U = at the critical dimension, implying no ultraviolet 
renormalization of A. This lack of renormalization holds to all orders in perturbation 
theory and follows from the non-anomalous character of the Galilei symmetry of the 
n — 1 model mentioned above. For an n = 1 driven diffusive system which is classically 



Galilei invariant but has a quantum mechanical anomaly, see ref. |L4] . 

To characterize the critical behaviour of the model for d < d c , we solve the renormal- 
ization group equations in the neighbourhood of all IR attractive fixed points 
Eqs. (H)-(||) form a system of differential equations describing the dependence of u, v and 
£ on /i. Its solution near an IR attractive fixed point gives, upon substitution in eqs. @- 
(tLOf), the /^-dependence of the anisotropy coefficients a and v near the fixed point. Indeed, 
in the neighbourhood of a fixed point the solutions to eqs. (|9|)-(|l0D read a oc // CT * and 
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v oc /A* , with 

C =T\ C =^±M 

1,0-* — J \* sv* — 1 -|- 2£ 

To have a signal for criticality, at least one of the two exponents C, a * and ( u * must be 
negative, since a or v will then diverge as // — > 0. To find the IR attractive fixed point 
with diverging length scales of the model for d < d c , we therefore proceed in two steps: 

Step 1. We first find all fixed points by solving the equations (3 U = f3 v = = and retain 
those which are IR attractive. We remind that for a fixed point to be IR attractive, the 
matrix of derivatives 



fiuu 


fiuv 


Put, 




fivv 


five, 




Ptv 





Pxy dy 



11) 



that results from linearizing eqs. ©-© about the fixed point must have eigenvalues with 
strictly positive real parts. The case — ► oo must be treated separately, since then there 
is no renormalization of £, hence no equation (||). In this case, fixed points are solutions 
of f3 u = (3 V = and the matrix of derivatives of the beta functions is the upper-left 2x2 
submatrix of (PI). 

Step 2. Next we calculate C, a * and for all the IR attractive fixed points found in step 1 
and keep those values which are negative. For finite, the condition (3% = and eq. (§) 
imply that the two critical exponents take the same value:^ Co-* = (u*- For — > oo, 
however, the condition (3% = is no longer required and Co-* and C^* m ay take different 
values. 

Although n has been assumed to take integer values in the range 1 to 4, the results 
can be analytically continued in n. In fact, analytic continuation gives well-defined beta 
functions for exotic values of n, like for example n —>■ or n —>■ —2. In what follows we 
take n as a real variable. 

In step 1 of this program, one must keep in mind that the values and v * of the 
couplings u and v at a fixed point are not physically meaningful. Hence redefining u = \ u 
and v = ipv, with \ and ip functions of £, does not change the physics. Such a redefinition 
may be used when one or both of w* and v* go to infinity to move the fixed point to m* and 

finite. The renormalization group equations for u and v are straightforward to obtain 
from eqs.([|)-([5p and involve the functions \ and ip. The functions \ and ip must be found 
such that the IR attractive fixed points of the new renormalization group equations have 

and finite. This introduces two new unknowns, x and ip, which, together with the 
complexity of the functions J 7 , Q, TC and C, makes it not feasible to analytically carry the 

3 This also follows from the definitions of v, £ a and Cn given in eqs. (||), (g) and (px|). 
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program above while keeping £ and n arbitrary. We have thus restricted ourselves to some 
values of n and £ that we find of interest. In particular, we have treated the cases 



1 n 



-2, 0, 1, 2, 3, 4 with arbitrary (w*, v t , £*), and 



(ii) = 0, (n — 2)/2n, oo with arbitrary (n, u* } v*), 

allowing in both instances for redefinitions of u and v as described above. The values 
= and = oo considered in case (ii) correspond respectively to 7i a d 2 s a and 7r a c^c^s^ 
being the dominant kinetic anisotropic term in the classical action S n . The idea of looking 
for critical points at = (n — 2)/2n stems from the simple relation between f vv and g vv 
for arbitrary n, namely 2n^g vv = {n — 2)f vv . Because of this relation, the equation (3^ = 
simplifies considerably if = (n — 2)/2n, thus making possible to find IR attractive fixed 
points with finite Another motivation to consider = (n — 2)/2n is to investigate 
whether the passage from n ^ 1 to n = 1 is smooth, since the case n = 1 can be recovered 
from = —1/2 and a — > oo. 

The results that we have obtained for (i) and (ii) are as follows (we provide the values 
of the physical quantities n, and Co-*, which characterize the universality classes). For 



n 






Co-* 




1/6 


-0.1257e 


-0.1257e 


3 


oo 


positive 


-0.2848e 




oo 


-0.1068e 


-0.4726e 




1/4 


-158e/629 


-158e/629 


4 


oo 


positive 


-0.1324e 




oo 


-0.3731e 


-0.3526e 



Table 1: Universality classes for n = 3,4. 

case (i) there are no IR attractive fixed points with diverging length scales if n = —2, 0, 2. 
If n = 1, the results of ref. |J are recovered. And if n = 3,4, we obtain the universality 
classes in Table 1. For case (ii) we obtain the universality classes: 

• £* = (n — 2)/2n, with n > 2 and critical exponents 



Co-* Cv 



Vv 



{n + 2) f vv + Ah 



(12) 



£=(n-2)/2n 



Note that the critical exponent is a complicated but well-defined function of n. 
• £* = oo, with n = n Q = (1 + y / 33)/4 and critical exponents 



Czv* o Co-* ^ 



(13) 
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= oo, with n > 2 and critical exponents 



± 1 , 2(n + l)X ± + nb 3 



where 



C* 2 nfe^ _ + hh _ hh (W) 

c ±_ 1 , u X ± [n(3n + 4)(n-2)6 1 -86 2 ]-8& 3 , . 

^* 16 1 J X± (6x6s -6264) +6166-6364 ' 1 ^ 



61 = — (3n 2 + 7n - 2) 64 = — (n 2 - 3n - 6) 

n — 2 n — 2 

b 2 = - - n (3n + 1) (n 2 + 2n - 12) 65 = - n (3n 2 - 3n - 10) 
8 2 

63 = - (2n 2 - n - 4) 6 6 = - - n 2 (n + 2) (2n 2 - n - 4) 



and 



1 



X± = — (-6 2 ± v /6 2 -46 1 6 3 ) . 

The exponents and are negative for all 2 < n < 4, whereas is negative only for 
2 < n < 2.0060. It is worth noting that the universality classes for n = —2,0,2,3,4 are 
particular cases of eqs. (|T2] ) -([T5| ) . 

The anomalous dimensions and found above are related to power laws in physical 
correlation functions and therefore should be directly measurable. In accordance with the 
continuity equation (H) and owing to the vector character of s a , physical correlators are 
expectation values of products of tensor currents 

JaP = ~0 d a Sf3 - 2(X£ d/sSa + ^ 5 a/3 s 2 + k (ns a sp - 5 a/3 s 2 ) . 

Since J Q p contains linear and quadratic terms in s a , the computation of physical correlators 
involves renormalization of composite operators. For n = 1, due to the non-renormalization 
of the coupling constant, renormalization of composite operators can be fully eliminated 
by considering correlators quadratic in the current |4j]. For n 7^ 1, since the coupling 
constants do renormalize, correlators cubic in J a p must also be considered. The analysis 
of physical correlators and their relation to the critical exponents ^* and C, a * is in progress 
and will be presented elsewhere. 

We conclude with a few comments: 

Comment 1. The anisotropy dimension n at which the critical exponents and 
occur is in general non-integer. It is particularly noticeable that for -u 2 = 0, v 2 finite and 
= 00, the only IR attractive fixed point with diverging length scales is that given in 
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eq. dl3|) and occurs at n = (1 + \/33)/4. With the caution due to the one-loop nature of 
our results, this suggests a fractal behaviour of the model. The similarity in some regards 
of our model with directed percolation may provide some insight in understanding these 
points. 

Comment 2. The fixed points associated with the critical exponents in eqs. (|12|)-(|l"3|) 
have either u* = or = 0, while the fixed points associated with the critical exponents 
in eq. (|14l)-(|T5"l) have u*, v* ^ 0. We have not been able to find IR attractive fixed points 
with all three parameters w*, and £* finite and non-zero. Our results show that if such 
a point exists to first order in e, it has fractal dimension. 

Comment 3. Eqs. (p^)-(pT5|) give for n > 2 two sets of critical exponents at £* = oo 
which are associated with two different attractive IR fixed points. This is so since ul 
and vl are two-valued functions of £ for £ — > oo and m*,^* ^ 0. Depending on the 
configuration of all fixed points, the system will flow in parameter space to one or another 
IR fixed point and will reach one or another set of critical exponents. 

Comment 4- If we take A = in S n , we find universality classes at £* = oo for n > uq, 
with critical exponents ( u * = and ( aif = — 2e/(n + 2). On the other hand, if we set 
k = 0, we obtain universality classes at £* = oo for n > 2, with critical exponents = 
and £o-* = — e/6. Note however that, despite the simplicity and beauty of these results, 
for A = or k = perturbation theory may generate at higher loops divergent relevant 
operators which are not in the classical action S n . For this reason, we have considered the 
more general case of two coupling constants. 
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